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Magnetite, Fe304, displays a highly complex low temperature crystal structure that may be 
charge and orbitally ordered. Many of the recent experimental claims of such ordering rely on 
resonant soft x-ray diffraction at the oxygen K and iron L edges. We have re-examined this system 
and undertaken soft x-ray diffraction experiments on a high-quality single crystal. Contrary to 
previous claims in the literature, we show that the intensity observed at the Bragg forbidden (00|)c 
reflection can be explained purely in terms of the low-temperature structural displacements around 
the resonant atoms. This does not necessarily mean that magnetite is not charge or orbitally ordered, 
but rather that the present sensitivity of resonant soft x-ray experiments does not allow conclusive 
demonstration of such ordering. 

PACS numbers: 



In many transition-metal oxides, the spatial localiza- 
tion of electrons on certain sites, so-called charge order- 
ing, has been used to explain some of their more in- 
triguing ground-state properties. For example, charge 
ordering has been invoked to describe phase transi- 
tions in some magnetoresistive manganites,'^' and the 
dynamic fluctuations of charge-ordered striped have 
been proposed as a mechanism of high temperature 
superconductivitjEl. Magnetite, Fe304, was the first ma- 
terial in which such a charge ordering transition was 
proposed, in connection with the metal-insulator tran- 
sition discovered by Verwey'*, and it has long been in- 
terpreted as the classic example of mixed-valence com- 
pound with formula unit Fe^+[Fe^+Fe^+]04 [Refs. I4|5j . 
In this interpretation, magnetite, which at room tem- 
perature crystallises into the cubic inverted spinel struc- 
ture AB2O4, with space group Fd3m, has formally Fe'^+ 
ions at the tetrahedral A sites and formally Fe^+ and 
Fe^"*" ions at the octahedral B sites. Unfortunately, this 
simple picture is deceptive as the crystal structure be- 
low Ty is extraordinarily complicated: the most recently 
reported structur^ consists of a complex P2/c mono- 
clinic cell containing 56 atoms in which the A and B Fe 
ions are split between two and six inequivalent sites re- 
spectively. Based on the P2/c structure, LSDA-fU band 
structure calculations have reported both charge order- 
ing (CO) (0.16 electrons) and an associate d 1 20 orbital 
ordering (00) on the octahedral sublattice^ET Seem- 
ingly arguing against this, however, are the results of 
resonant x-ray scattering (RXS) experiments at the iron 
K-edge, have been interpreted a s provi ding evidence ei- 
ther against any charge orderinJ^HnUI] or in favor of a 



0.12 electrons charge disproportionatiorP^ between the 
formally Fe^+ and Fe'^+ sites than predicted in Refs.'^. 
Very recently two further papers have appeared in which 
soft-x-ray diffraction measurements at the O K-edg^^^ 
and the Fe L-edge^^ were interpreted as providing evi- 
dence for both charge and orbital ordering. We not that 
in the monoclinic cell the iron B sites are no longer equiv- 
alent by symmetry and there is therefore no requirement 
that they have the same charge density surrounding the 
atomic site. It is thus most likely that they will not have 
the same charge density. The question is therefore: what 
is the smallest charge difference, about which one would 
reasonably claimed that the material is charge ordered? 

Resonant x-ray scattering occurs when a photon ex- 
cites a core electron into an excited state and is sub- 
sequently re-emitted when the electron and core hole 
combinelSl. On resonance the x-ray scattering amplitude 
is anisotropic and is sensitive to the anisotropic charge 
distribution of the resonating ion. The anisotropic charge 
distribution can be intrinsic to the scattering ion due 
to orbital occupation or can be intrinsic to the lattice 
as in the case of Templeton-Templeton scatterin^i^EZIIll 
The characteristic of Templeton-Templeton scattering is 
that a reflection which is Bragg-forbidden because of a 
compound symmetry operation, such as a ghde plane or 
screw axis, becomes allowed when the incident photon 
energy is tuned to a resonance. On resonance the x-rays 
are sensitive to the quadrupolar term in the charge dis- 
tribution, of the resonating atom, Q, and the difference 
between the two electric quadrupole moments, related 
by the symmetry operation, sum to zero and a resonant 
peak is observed arising from the crystal structure. 
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We have chosen to revisit magnetite and report here on 
resonant soft x-ray experiments that confirm the resonant 
enhancement of the (OOi)c reflection at both the oxygen 
K- and iron L3 edges. However, we have carried out a 
careful analysis of this superlattice reflection focussing 
on a detailed investigation of the effects of the distorted 
crystal structure below the Verwey transition, without in- 
voking any charge or orbital ordering. We find that we 
can model our data well by considering only Templeton- 
Templeton scattering arising due to the structural distor- 
tions below the Verwey transition, without the need to 
resort to charge order or orbital order. This is contrary 
to the claims of Refs. I13I14I 

The experiments were conducted on high-quality syn- 
thetic magnetite crystals prepared in an arc-image fur- 
nace using the floating-zone technique. The purity of 
the sample was verified by heat capacity measurements, 
which gave a maximum heat capacity value of 120.6 K 
and an entropy change of 5.77 JK~^ at the Verwey tran- 
sition. These results give for Fe3_504 a value S = 0.0002 
showing that the stochiometry of our sample is very close 
to the ideal case. Soft x-ray diffraction experiments were 
conducted on the ID08 beamline at the ESRF, Greno- 
ble, France. In what follows we index the sample in the 
approximate low temperature Pmca orthorhombic struc- 
ture (No. 57) with lattice parameters a — 5.944 A, 
h = 5.925 A and c = 16.775 A^l. This structure is re- 
lated to the P2/c (No. 13) structure by only a slight 
monoclinic distortion (/3 = 90.2363°). In this orthorhom- 
bic setting, the cubic (00|)c reflection becomes the or- 
thorhombic (001)0- The sample of Fe304 was cut with a 
[001]o surface normal and polished with 0.1 \an diamond 
paste. It was then mounted on a SmCo magnet, provid- 
ing a field at the sample surface of « 0.3 T, parallel to the 
surface normal. This field defines a unique c-axis so that 
on cooling through the transition the number of crystal- 
lographic domains are minimized. The sample was then 
to a base temperature of 30 K and a resonant signal was 
observed at the (OOl)o position in reciprocal space in the 
vicinity of the iron and oxygen K-edges. 

Figure [T] shows an incident photon energy scan at con- 
stant wavevector of (OOl)o as the is tuned in the vicin- 
ity of the oxygen K-edge. The scattering is observed to 
peak at 529.1 eV, about 10 eV below the main oxygen 
K-edge. The insert shows a scan along the [001] o di- 
rection through the (OOl)o reflection at Ei = 529.1 eV, 
with a flt to a Lorentzian squared lineshape. The cor- 
relation length obtained from this flt is > 3000 A. This 
represents a lower bound on the penetration depth of the 
x-rays and thus this value indicates that the resonant sig- 
nal is not surface sensitive even at the maximum of the 
resonance. The (OOl)o reflection was also visible in the 
vicinity of the iron L3 edge (Fig. |2| . The bottom panel 
of Fig. |2] shows the incident photon energy dependence of 
the integrated intensity of the (OOl)o reflection. The ex- 
perimental signal is only visible at the L3 threshold, with 
a maximum at an energy of 706.5 eV and is found to be 
suppressed above 708 eV. Such behavior arises from the 
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FIG. 1: (Color online) The incident photon-energy- 
dependence of the (OOl)o reflection in Fe304 close to the oxy- 
gen K-edge (red circles). The solid black line represents ab- 
initio calculations of the scattered intensity assuming struc- 
tural distortions and no charge or orbital order. In the insert, 
a scan along the [001]o direction through the (OOl)o reflection 
is shown (red circles); the solid line is a fit to a Lorentzian 
squared lineshape. 



very large self-absorption caused by the strong Fe L3 res- 
onance, leading to a total loss in the observed signal. The 
width of the diffraction peak as a function of energy as 
shown in the top panel of Fig. [2]conflrms this. The peak 
width is broader than that found at the oxygen K-edge 
and tracks the calculated absorption (dashed line) indi- 
cating that the change in width arrises from the increased 
absorption and consequently reduced penetration depth. 
Finally, Fig. [3] shows the temperature dependence of the 
integrated intensity of the (OOl)o reflection measured at 
both the iron L3 and oxygen K edges. The data were col- 
lected by performing rocking scans of the sample angle, 
6, at each temperature. The signal at both edges was 
found to be virtually constant up until a temperature of 
w 125 K above which no intensity is observed. 

We now turn to our resonant scattering simulation. 
We have used the FDMNES prograrrPSl in the multiple 
scattering mode. *** The results of these sumulations 
are shown in Fig. [l] and Fig. [2] for the oxygen K and iron 
L edges respectively. In order to calibrate to the experi- 
mentally obtained data with the FDMNES simulations, 
the calculated absorption was compared with the sample 
absorption measured by total electron yield at the oxygen 
K-edge. Our simulation reproduces well the main experi- 
mental features, including the energy gap of about 10 eV 
between the RXS signal and the main oxygen absorption 
edge, as well as the energy width of the peak. 

In this speciflc case, there are eight in-equivalent oxy- 
gen sites of Ad Wyckoff symmetry in the Pmca space 
group. Considering only the oxygen atoms which dom- 
inate at this energy, the structure factor of the (OOl)o 
reflection is: 

•^(001)0= X! 2/,(l - mj,)cos(27rWj), (1) 

] = 1...8 
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where fj is the atomic scattering ampHtude, {j = 1 . . .8 
labels the inequivalent sites), Wj is the fractional coor- 
dinate of the jth oxygen atom in the c direction and 
rhy is the mirror plane in the b direction of the Pmca 
setting. Using the local mirror symmetry ihx of the 4(i 
sites, we find that fj oc Q^y^, the electric quadrupole ma- 
trix element. In the monoclinic P2/c setting, the fa^ 
symmetry is lost and a further contribution fj cx Qiy 
appears. This can be shown to be negligible, since it is 
proportional to the small angular distortion, (3 ^ 90°, 
from the orthorhombic Pmca structure. Upon evaluat- 
ing the structure factor S we can conclude that almost 
all the scattered intensity comes from the sum of the 
quadrupoles Qyz at oxygen sites Oi and O2 only. That 



is S, 



(ooi)„ 



Qyz + Qvz high-temperature phase 



iyz 



Qy^ = —Qyz, due to the C2Z screw axis of the high tem- 
perature Fd2>m space group, and therefore their sum is 
zero. Below the Verwey transition, the unequal atomic 
displacements of the Oi and O2 oxygen sites from their 
high temperature positions, gives a finite signal even if 
Qy^ = —Qy^. However, this signal is tiny because of the 
very small displacement, and has an expected amplitude 
of K, 10^^ X IQ^^^I- contrast, a much bigger amplitude 
might be expected if the surrounding iron tetrahedra are 
distorted making Qy^ ^ ~Qyz ■ 

To investigate which of these contributes to be the 
most significant we have performed several numerical cal- 
culations in which the tetrahedral or octahedrally coordi- 
nated iron atoms, and/or the oxygens were in turn placed 
in their high-temperature positions, with the rest of the 
cluster held in their low-temperature positions. By this 
method, we found that the main contribution to the sig- 
nal comes from the O2 position (~ 70% of the total): in 
fact Q'^^ does not change much when the iron sites move 
from the high-temperature to the low-temperature posi- 
tions, while Qy^ varies by about 100 %. This dominant 
change is due in particular to the displacement of the 
octahedral iron atoms surrounding the O2 cite. The iron 
sites belong to the Fe i?3 sites, in the notation of Ref. 
and undergo the strongest distortion when passing from 
the high-temperature to the low-temperature phase: the 
FeB3-02 distance changes from 2.06 A to 1.96 A, about a 
5% contraction. Therefore the O K-edge signal is mainly 
determined by the hybridization of 2p oxygen orbitals at 
O2 sites with 3d iron orbitals belonging to octahedral Fe 
B'd sites. 

The fact that we can explain the signal through these 
atomic displacements contradicts the interpretation of 
the signal at (OOl)o in Ref. TS", that explicitly excluded a 
structural origin to the (OOl)o reflection at the oxygen K- 
edge. In particular, the arguments made by the authors 
of Ref. [131 for the assignment of their signal to charge 
and orbital order due to its polarization dependence, are 
reproduced in our calculations based solely on structural 
distortions. This makes clear that invoking charge or or- 
bital ordering to explain the detection of this superlattice 
reflection is unnecessary, and potentially misleading. 

This same procedure was then used to evaluate the res- 
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FIG. 2: (Color online) Top panel: Measured total electron 
yield (x-ray absorption spectrum) at 27 K (red dashed line) 
and calculated absorption spectrum (solid black line) through 
the Fe L2,3 edges. Bottom panel: incident photon energy de- 
pendence of the integrated intensity of the (OOl)o reflection 
close to the Fe L-edges (red circles). The solid black line 
shows ab-imtio calculations of the resonant scattering (see 
text). The inset represents a scan along the [001] o direction 
at an energy corresponding to the maximum in the signal. 
The variation of the longitudinal width of the (OOl)o reflection 
with incident energy is included in the top panel for compar- 
ison with the experimental and calculated absorption (blue 
diamonds) . 



onant signal of the (OOl)o reflection at Fe L2.3d-edges, as 
shown in Fig. [2] The calibration between our experimen- 
tal data and the FDMNES simulation was set by com- 
parison of the calculated absorption with the absorption 
as measured by total electron yield. 

We can repeat the same analysis as performed at the 
oxygen K-edge for the data at the iron L2,3 edges. The 
results of the FDMNES simulation are shown in Figure |2] 
In the top panel a comparison between the measured and 
calculated absorption are shown, while the bottom panel 
shows the comparison between the integrated intensity 
of the (001)0 reflection as a function of incident photon 
energy and our simulation. 

In the Pmca setting there are 6 inequivalent groups 
of iron atoms, with two groups of tetrahedral iron-sites 
{Al and A2, following the notation of Ref. [20 ), and 4 
groups of octahedral iron-sites {Bl, B2, B3, and 54), 
each group containing 4 iron atoms. Al, A2, B3, and 
BA sites have a local m^^-symmetry {Ad Wyckoff site), 
so that the same considerations discussed above for the 
structure factor at the oxygen K edge are still valid. The 
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FIG. 3: (Color online) Temperature dependence of the (001) 
reflection of Fe304 at the O K-edge (red circles) and the iron 
L3 edge (blue diamonds) . 



only signal that can be measured for these ions is due to 
the this time projected on the corresponding iron- 
sites j. The B2 site has a C2y local symmetry (4c Wyck- 
off site) and the two groups of two ions that contribute 
in antiphase at the (OOl)o are related by inversion sym- 
metry, so their total contribution equal to zero. Finally 
the Bl sites, with local inversion symmetry (46 Wyck- 
ofF site), also contribute with the quadrupole component 
Qijz- When we numerically compare the separate con- 
tributions of the Al, A2, Bl, B3 and B4 sites, we find 
that the intensity from the Bl site is smaller by a factor 
of 500 relative to the contribution from the 34 site. We 
can conclude therefore, that the dominant contribution 
of the intense low-energy peak seen in the calculations 
at 708.5 eV is from the iron atoms sitting at the B4 
sites. The tetrahedral A sites and the octahedral B3 site 
contribute mainly to the smaller shoulder at higher en- 
ergy, and the fact that they have opposite amplitudes 
gives rise to the local minimum at 710 eV. Therefore all 
the experimentally detected signal comes from the B4 
sites (nominally Fe^+). This is in keeping with some pre- 
vious results obtained with rather different approaches 



which concluded that the low-energy part of L3 spec- 
trum (around 707 eV) is mainly determined by the t2g 
states of the nominally divalent iron ion j^ l ^^ l However 
there is an important difference between the work pre- 
sented here and atomic-multiplet-based calculations such 
as those reported in Refs. I14|21l By focussing on the or- 
bital occupancy in a ionic model, these latter authors 
have neglected the important structural differences that 
exist between sites like Bl and Bi (that are otherwise 
equivalent where formal charge is concerned, both being 
formally Fe^"*"). The fact that we calculated their relative 
contribution to the total signal in the ratio 1:500 proves 
such a difference. 

The similar temperature dependence shown in Fig. [3] 
can now be explained as a natural consequence of the 
signals at the oxygen K-edge and iron L-edges both mea- 
suring the same order parameter. As discussed above we 
argue that this order parameter is structural distortions 
associated with the structural phase transition from the 
Fd3>m high temperature structure to the low tempera- 
ture P2/ c structure and not that of any charge or orbital 
order. 

In conclusion, we have shown that the (OOl)o reflection 
of Fe304 is sensitive to the local displacements around 
the resonant ion at both the oxygen K and iron L edges. 
The electronic anisotropy arising from the crystal distor- 
tions are sufficient to explain the origin of the scattered 
signals, and any invocation of charge ordering and/or or- 
bital ordering is not necessary to reproduce the data. At 
the oxygen K-edge, the signal is determined by the hy- 
bridization of O 2p orbitals of the four O2 atoms with 
the B3 Fe 3d orbitals (nominally Fe'^^). At the iron L3 
edge, in contrast, the resonant x-ray scattering is mainly 
sensitive to the contribution of Fe 3d orbitals from BA 
sites (nominally Fe^+). 
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